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We demonstrate precise control over the zero-phonon optical transition energies of individual nitrogen- 
vacancy (NV) centers in diamond by applying multiaxis electric fields, via the dc Stark effect. The Stark shifts 

display surprising asymmetries that we attribute to an enhancement and rectification of the local electric field 
by photoionized charge traps in the diamond. Using this effect, we tune the excited-state orbitals of strained NV 
centers to degeneracy and vary the resulting degenerate optical transition frequency by >10 GHz, a scale com- 
parable to the inhomogeneous frequency distribution. This technique will facilitate the integration of N V-center 
spins within photonic networks. 

PACS numbers: 71.55.Cn, 71.70.Ej, 78.56.-a, 81.05.ug 



Nitrogen-vacancy (NV) centers in diamond are promising 
solid-state qubits for emerging quantum technologies, due to 
their long spin coherence times [1] and fast manipulation rates 
[2], together with a level structure that allows for straight- 
forward optical initialization and readout of the electronic 
spin state [3]. Furthermore, in high-quality single-crystal di- 
amond at temperatures below «25 K [4], sharp zero-phonon- 
line (ZPL) optical transitions facilitate the coherent coupling 
between NV-center spins and photons [5, 6]. The integration 
of NV centers within photonic structures [7] to route single 
photons and enhance the spin-photon interaction could there- 
fore lead to scalable applications for quantum iirformation 
processing and secure communication [8]. 

As solid-state "trapped atoms," NV centers are sensitive to 
their local environment. While this sensitivity has enabled 
nanoscale magnetic [9] and electric [10] metrology, it also ex- 
poses individual NV centers to sample inhomogeneities, lead- 
ing to a distribution of ZPL frequencies within a diamond 
[11]. The ability to tune these frequencies is crucial for pho- 
tonic applications, for instance to utilize the selection rules at 
the symmetry point for spin-photon entanglement [6] or 
to coherently couple distant NV centers to indistinguishable 
photons. Through the dc Stark effect, applied electric fields 
perturb both the ground-state spin [10, 12] and excited-state 
orbitals [13], providing the means to control the optical tran- 
sitions. 

Here we use micron-scale devices to manipulate electric 
fields in three dimensions, to compensate the intrinsic local 
strain and electrostatic fields of individual NV centers and 
achieve full control of the orbital Hamiltonian. Furthermore, 
by analyzing the Stark shifts as a function of applied volt- 
ages, we infer a surprising amplification and rectification of 
the local electric field, consistent with electrostatic contribu- 
tions from photoionized charge traps within the diamond host. 
By harnessing this reproducible effect, we can tune the NV- 
center Hamiltonian to arbitrary points across a range compa- 
rable to the inhomogeneous ZPL distribution. 

The electronic structure of the negatively charged NV cen- 
ter is determined by symmetry, through its point group C^y 



[11, 14]. The spin-triplet ground (GS, symmetry ^A2) and ex- 
cited states (ES, symmetry ^E) are connected by ZPL tran- 
sitions around 637.2 nm (1.946 eV). Our experiments are 
performed at zero magnetic field, where the spin-triplet basis 
states are {l^^) , l^j,) , l^^)}. A 532 nm (2.3 eV) "repump" 
beam pulsed at ~300 kHz in a confocal geometry maintains 
a spin-polarized population in \Sz). Between repump cycles, 
we count photoluminescence excitation (PLE) photons emit- 
ted by the NV center into the redshifted phonon sideband [see 
Fig. 1(a)] after absorption from a narrow-line red laser tun- 
able across the ZPL transitions. As we scan the red laser fre- 
quency, we typically measure two peaks in the PLE spectrum 
as shown in Fig. 1(b); these correspond to spin-conserving 
transitions from the GS orbital singlet \A2,Sz) to the two ES 
orbital eigenstates {|Ei, S^) , |E2, Sz)}. In a crystal environ- 
ment with perfect C^y symmetry and zero electric and mag- 
netic fields, these ES orbital states would be degenerate, but 
the symmetry is generally broken by local crystal strain and 
by nonuniform electrostatic charge distributions that generate 
local electric fields. 

The dc Stark perturbation to the Hamiltonian, Hstark = 
fi F, describes the interaction between the local electric 
field F and the electric dipole operator fi. For fixed stress, 
the strain perturbation can be cast into the same form by iso- 
lating components which transform as the irreducible repre- 
sentations of C^y [14]. The combined perturbation has the 
form VaiOai + Ve^Oe, + Ve^Oe^, where Or„ is an orbital 
operator transforming as the basis state ir^) and 

f Va, = Sa, - l^\\F, 

[ Ve^ = 5e„ - ii±Fy 

are the symmetrized field strengths, in terms of fixed strain 
components ^r^, projections of the local electric field Fj, 
and the reduced matrix elements of the electric dipole oper- 
ator {/U||,/Uj_}. We choose orbital basis states {|Ea;) , |Ey)} 
for the ES which transform like vectors {x,y} in the NV- 
center coordinate system [15 1, and we ignore the small 
(pslOO MHz) spin-spin coupling between ES spin states \Sz) 
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FIG. 1 . (a) Simplified energy-level diagram (not to scale) showing 
only the l^z) levels of the NV-center ground and excited states, with 
resonant excitation (red arrows) and redshifted emission (gray arrow) 
marked, (b) PLE spectrum (points) with no applied bias, marked 
by dashed line in (d), with a two-Lorentzian fit (solid curve), (c) 
Micrograph and photoluminescence image (center) of device A, with 
electrical connections marked, (d) PLE spectra of the 6-|J.m-deep NV 
center marked by an arrow in (c) as a function of lateral bias applied 
symmetrically to the X (left panel) or Y (right panel) gate pairs, with 
Vdc = V. In all PLE spectra, the origin of the relative frequency 
axis is arbitrary. 

and [16]. By defining ^|A2,5,) = 0, the 

Hamiltonian in the {|E^, iS^) , |Ej,, S'^)} basis can be written 
as 

i7^(;^o + A.„^.)I+^(_\ i;:^"), (2) 

where hujo is the natural transition energy including fixed per- 
turbations of Ai symmetry, and A/i|| = — defined 
such that both A/xy and fj,± are positive. The transition energy 
eigenvalues take the form E± = hD ± ^hS, where 

/li? = fiwo + A/i||Fj, (3a) 

h5^V2(vl+Viy^^ (3b) 

are the longitudinal and transverse components due to fields 
of Ai and E symmetry, respectively. From Eq. (1), it is clear 
that a local electric field can cancel the transverse components 
of intrinsic strain to restore C^y symmetry to the system, and 
from Eq. (3a) we see that an electric field applied along the 
NV-center symmetry axis shifts the energy of both transitions 
by the same amount. 

We first investigate these effects using device A, shown in 
Fig. 1(c), consisting of four Ti-Pt-Au gates fabricated on the 



diamond surface. The sample is a 0.5-mni-thick single-crystal 
diamond grown by chemical vapor deposition with <5 ppb ni- 
trogen content (ElementSix), irradiated with 2 MeV electrons 
(1.2 X 10" cni-2) and then annealed at 800 °C to create NV 
centers. Measurements are performed in a continuous-flow 
cryostat operating at «20 K. Symmetric biases Vx and Vy 
applied as shown in Fig. 1(c) produce lateral electric fields 
Fx and Fy in the [110] and [110] crystal directions, respec- 
tively, while a common dc bias generates fields in the [001] 
out-of-plane (Z) direction. The sample {X, Y, Z) and NV- 
center {x,y, z) coordinate systems are uniquely related for a 
given NV-center projection from the (111) family [17]. 

Figure Id contains a series of PLE spectra showing the opti- 
cal resonances of the NV center marked in Fig. 1(c), as a func- 
tion of separate biases Vx and Vy with Vdc = V. The lateral 
biases are applied as symmetrically pulsed square waves at 
1 kHz, with PLE photons binned according to polarity. While 
the response to static lateral bias is qualitatively similar, this 
technique separates slow photoinduced charging effects from 
the dielectric response as discussed below. Two features are 
evident in the data: first, we observe an unexpected "kink" 
at zero bias, and second, the resonances cross at Vx « 7 V, 
demonstrating that we can indeed restore Ca^ symmetry to 
the system. We explore both of these features below with ad- 
ditional experiments. 

The kink at zero bias reflects an asymmetry in the local 
electric field vector as a function of polarity, i.e., F{+V) ^ 
—F{~V). We argue that this asymmetry results from the pho- 
toionization of charge traps in the diamond host. Even in 
high-quality single-crystal synthetic diamonds, deep defects 
such as vacancy complexes and substitutional nitrogen have 
important effects on the material's electronic properties [18]. 
In particular, substitutional nitrogen atoms form donor levels 
sal.7-2.2eV below the conduction band edge [19], and the 
timescale for charge transport through these levels is very long 
(hours) even in nitrogen-rich diamond at room temperature 
[20]. These traps are easily ionized by the 532 nm repump 
beam (salOO M-W) which is 4-5 orders of magnitude stronger 
than the red laser («1 nW). When voltages are applied, this 
leads to a long-lived nonequilibrium charge distribution in the 
illuminated volume of the sample, which can either amplify 
or screen the local electric field. 

As a simplified one-dimensional demonstration, we present 
in Fig. 2(a) the Stark-shift response of an NV center 13 \xm 
below the surface of a second diamond sample irradiated and 
annealed under similar conditions, but patterned with a global 
top gate of the transparent conductor indium-tin-oxide (device 
B). As the top-gate bias is stepped in a loop over Ral60 min, 
we observe hysteresis in the response characteristic of wl h 
charge-relaxation timescales. Furthermore, by comparing the 
magnitudes of the Stark shifts due to biases applied later- 
ally across an 8 ixm gap [Fig. 1(d)] and vertically across the 
0.5 mm sample thickness [Fig. 2(a)], we find that, when the 
top-gate bias is negative, the local electric field below the top 
gate appears to be amplified by roughly an order of magnitude 
over dielectric predictions [17]; conversely, the field appears 
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FIG. 2. (a) Stark-shift hysteresis loop for an NV center 13 |J.m be- 
low a transparent top gate, as a function of top-gate voltage. Points 
mark the transition frequencies from a two-Lorentzian fit to a PLE 
spectrum at the corresponding voltage and color-coded arrows indi- 
cate the sweep direction, (b) Schematic of the local electric fields in 
a lateral geometry. Photoionized charge traps in the illuminated vol- 
ume contribute a rectified field R predominantly in the +Z direction, 
which adds to the dielectric field E to shift the direction of the local 
field F. (c) dc Stark components v and 5 (points) extracted from fits 
of the PLE spectra in Fig. 1(d), with a combined fit according to the 
model described in the text (solid curves), (d) dc Stark components 
(points) and fits (solid curves) measured both with (green squares) 
and without (orange circles) the 532 nm repump excitation. Inset: 
Micrograph and photoluminescence image of device C, with electri- 
cal connections marked. The 7-|j.iii-deep NV center measured in (d) 
is circled. In all cases, marker sizes slightly exceed measurement 
uncertainties. 

to be completely screened above a threshold bias, where the 
response is flat. The response of an NV center in device A 
to variations of Vdc is qualitatively similar [17], and both are 
consistent with a picture in which positive charges in the illu- 
minated volume below the NV center rectify the Z component 
of the electric field. 

We can incorporate these charging effects into a phe- 
nomenological model capturing the essential features of our 
observations. As depicted in Fig. 2(b), the local electric field 
for an NV center between two surface gates is composed of 
a dielectric component roughly parallel to the sample surface 
and a rectified component due to photoionized charge that is 
mainly out of plane. We model this field as 

F = /?Fv + /3|y||, (4) 

where /3 accounts for geometric and dielectric factors that pre- 
dict a local electric field in the direction v in response to 



an applied voltage V , and ^ is a dimensionless vector giv- 
ing the relative strength and direction of the rectified field, 
assumed to scale linearly wifli \y\. Because of the long 
charging time scale, the rectified field /3|V^|^ does not change 
when we switch the bias polarity on millisecond time scales 
while the dielectric component changes sign, produc- 
ing a polarity-asymmetric response. The assumption of a lin- 
ear relationship between the rectified field strength and is 
motivated by the empirical observation that P, proportional to 
Fz, varies linearly with applied bias in all our measurements. 
This amounts to an approximation that the spatial distribution 
of photoionized charge remains fixed, while the charge den- 
sity varies linearly with \y\. 

Figure 2(c) shows the mean (p) and difference {S) of the 
transition frequencies extracted from fits to the PLE spectra in 
Fig. 1(d). The NV-center symmetry axis ([111] in this case) 
is uniquely determined by the sign of v in response to elec- 
tric fields in different directions. By substituting Eq. (4) into 
Eqs. (3) and applying the appropriate coordinate transforma- 
tion, we obtain a model that quantitatively agrees with our 
observations [17], as shown by the fits to the data in Fig. 2(c). 
Given that this is only a simplified phenomenological descrip- 
tion of a complicated three-dimensional system, it matches 
our observations surprisingly well. 

Finally, we present a control experiment in which we miti- 
gate effects due to the 532 nm repump cycle. Occasional re- 
pump pulses are still required to compensate for photoioniza- 
tion of the NV^ charge state due to sequential two-photon ab- 
sorption, but with weak (<1 nW) resonant light, the required 
repump period can be increased to several seconds [4], allow- 
ing time to apply bias, record a complete PLE spectrum, and 
rezero the bias, all between repump pulses. Since weak spin- 
nonconserving optical transitions quickly polarize the NV- 
center spin away from resonance in the absence of the repump 
cycle, we mix the spin population by applying a microwave 
magnetic field resonant with the GS spin transition. We use 
another device for this purpose (device C), shown in Fig. 2(d), 
which is fabricated on the same diamond as our four-gate lat- 
eral device. It consists of a short-terminated waveguide to 
generate microwave fields and serve as ground, and a gate that 
when biased produces lateral electric fields across an 8 jam 
gap- 
Figure 2(d) shows V and (5 as a function of gate voltage 
for the NV center circled in the inset. Once again the bias 
polarity is switched at 1 kHz and the PLE photons are binned 
accordingly. A polarity asymmetry is clearly observed when 
the repump beam is present, particularly as a kink in v, and 
it is significantly reduced when the biases are applied in the 
absence of the repump cycle. Fits to the data using our model 
[17] are shown as solid curves, from which we find that |^| 
is reduced from 0.71 ±0.02 with the standard repump cycle to 
0.33±0.03 when the 532 mii beam is omitted. 

Based on this understanding, we can exploit the photoin- 
duced charge to obtain greatly enhanced tunability in our four- 
gate geometry (device A). Since the rectified field points pre- 
dominantly out of plane and has strength comparable to the 
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FIG. 3. Tuning diagram for the NV center marked in Fig. 1(c). 
Degeneracy is achieved at different frequencies by setting the ap- 
plied biases (Vdc, Vx, Vy) as marked in the upper panel; in the lower 
panel we show PLE spectra as a function of Vx around each of these 
points for fixed Vdc (Vy is also varied to keep the ratio Vx /Vy con- 
stant). We can shift a second NV center to degeneracy within the tun- 
ing range of the first, as shown by the PLE spectra outlined in blue 
(lower panel) and the corresponding blue line (upper panel) marking 
the degenerate frequency. 



dielectric component, we effectively obtain three-dimensional 
control of the local electric field vector. The application of a 
negative reference bias Vdc as shown in Fig. 1(c) increases 
the rectified component Fz independently of {Fx, Fy). As a 
demonstration of the flexibility of this technique, we present 
a tuning diagram in Fig. 3 in which we use {Vx, Vy) to scan 
through the C^v degeneracy point at different settings of Vdc. 
Each crossing occurs at a different frequency, with the corre- 
sponding bias point {Vx , Vy , Vdc) marked in the upper panel. 
Essentially, we are compensating the transverse components 
(S'e^ , <S'e^) of the intrinsic fields and tuning the longitudinal 
component of F to shift the frequency. Since the rectified 
field always points along +Z, we can only tune the frequency 
in one direction, but the effect is strong enough to produce 
a >10 GHz shift in the degenerate frequency with practical 
applied voltages. 

With this technique we can tune multiple NV centers to 
have the same degenerate transition frequency. The PLE spec- 
tra outlined in blue in Fig. 3 were obtained from a second 
NV center in the same device at Vdc = V and display 
C^y degeneracy at a frequency within the tuning range of the 
first. If these two NV centers were in separately controlled 
devices and tuned simultaneously to degeneracy at the same 
frequency, they would couple identically to indistinguishable 
photons. 



In conclusion, we have used electric fields to tune the 
ZPL transitions of individual NV centers in micron-scale de- 
vices size compatible with photonic structures. Through their 
dc Stark shifts, NV centers serve as nanoscale probes of 
their electrostatic environment, revealing strong signatures of 
charge accumulation due to photoionization of deep donor 
levels in the diamond. We have analyzed these effects with 
a phenomenological model and used the additional fields pro- 
vided by photoionization to obtain three-dimensional control 
of the local electric field, in order to tune both the overall en- 
ergy and orbital splitting of the excited-state Hamiltonian. In 
particular, we have demonstrated how to reach the C^v sym- 
metry point and then apply longitudinal perturbations to shift 
the degenerate photon energy. By coupling multiple NV cen- 
ters to indistinguishable photons with this technique, photonic 
networks could provide a quantum bus to coherently couple 
distant NV centers, and entanglement swapping protocols [8] 
could enable long-distance quantum key distribution. 
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METHODS AND SUPPLEMENTARY DATA 

Optical methods 

Measurements were performed using two home-built con- 
focal microscopes designed to excite single NV centers both 
resonantly at 637.2 nm and in the blue-shifted phonon side- 
band at 532 nm, while collecting photoluminescence photons 
from the red-shifted photoluminescence sideband between 
R:i650-830 nm. Details and a schematic are available else- 
where [1]. The 532 nm laser serves both as the excitation 
for standard photoluminescence spatial imaging and as a 're- 
pump' cycle to maintain the NV^ charge state which is pho- 
toionized by resonant excitation alone. 

Resonant excitation is provided by tunable diode lasers 
from the New Focus Velocity^'^ series, which are tuned by 
means of a piezoelectric actuator that makes fine adjustments 
to the laser cavity. An important technical aspect of these ex- 
periments is the calibration of the laser frequency response 
to variations in the voltage applied to the piezo, since this is 
typically nonlinear and hysteretic. In one of our microscope 
setups we have pre-calibrated the laser frequency response to 
particular piezo modulation sequences using a Mach-Zehnder 
interferometer. In the other setup, unpredictable mode hops 
of the laser make such calibrations unreliable, so we con- 
tinuously monitor the laser frequency interferometrically: a 
Fabry-Perot cavity with a free spectral range of 1.5 GHz pro- 
vides fine frequency resolution, while a Mach-Zehnder inter- 
ferometer with free spectral range ss7 GHz is used to detect 
mode hops. In both cases, our calibrations provide relative 
frequency determinations with salOO MHz uncertainty across 
the «90 GHz piezo modulation range. The resonant light is 
prepared with circular polarization, so that the (linear) opti- 
cal selection rules for 15*2) optical transitions do not cause the 
resonances to disappear in certain regimes. 

In both setups, we also have the ability to add confocal exci- 
tation at 405 nm (3.1 eV) to the optical path. Unsurprisingly, 
this higher-energy excitation significantly alters the charging 
dynamics of the system. Generally it causes the NV center 
optical transitions to become less stable and so it was not used 
during any of the measurements presented in the main text, but 
we did use it as a control to 'reset' the equilibrium charge dis- 
tribution with no applied bias between experiments. A short 
burst of 405 nm light applied in this configuration was gener- 
ally enough to restore the initial NV center transition frequen- 




FIG. SI. Stark shift hysteresis loop for the NV center marked in 
Fig. Ic of the main text as a function of Vdc (with Vx, Vy ~ 0). 
Points mark the optical resonance frequencies extracted from two- 
Lorentzian fits to PLE spectra at the corresponding voltage, and 
color-coded arrows mark the sweep direction. Errorbars represent 
the frequency uncertainties from the fits. 



cies (wlO s at optical power wlO ixW), while full relaxation in 
the presence of the 532 nm light alone required several hours. 



Charging behavior in device A 

The interpretation of the hysteresis measurements pre- 
sented in the main text (Fig. 2a) is simplified by the one- 
dimensional geometry of the top-gate device, but we mea- 
sured a qualitatively similar response in our four-gate lateral 
device (device A, see Fig. Ic in the main text) to variations of 
the reference voltage Vdc common to all four gates. Figure 
S 1 shows a hysteresis loop as a function of Vdc for the NV 
center marked in Fig. Ic, which is w7 |J.m below the diamond 
surface. Although the response is similar to the top-gate mea- 
surements we do notice a few differences: First, the hysteresis 
'direction' is reversed, in the sense that the Stark response 
lags Vdc in Fig. SI, while it leads the top-gate bias in Fig. 2a. 
Second, the 'threshold bias' is closer to zero in the lateral de- 
vice than in our top-gate measurements. These differences 
are probably due to a variety of factors, including (a) geomet- 
ric differences due to the gate geometry and the depths of the 
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respective NV centers, (b) variations in the diamond samples, 
(c) details of the diamond-metal interfaces for different met- 
als (Ti/Pt/Au on the lateral device and indium-tin-oxide as the 
top gate), and perhaps most importantly (d) the electrostatic 
boundary conditions at the back of the device - the top-gate 
device B is indium-bonded to ground on the backside, while 
device A is glued to the ground plane with an insulating adhe- 
sive. 

We note that polarity-dependent Stark shifts have been re- 
cently observed for Chromium centers in diamond [2], and 
suggest that a one-dimensional experiment with a transparent 
top gate could determine whether these devices display simi- 
lar charging behavior 



Photoconductivity 

When voltages are applied between surface gates in the 
presence of optical excitation, we expect the photoexcitation 
of defect levels to lead to a photocurrent near the surface. 
Fig. S2 shows the photoconductivity response of device A 
when the excitation beam is focused on the NV center marked 
in Fig. Ic of the main text, 6 iJ.m below the diamond surface. 
The top panel shows an example of a typical dc photoconduc- 
tivity response for this and similar devices we have measured, 
with almost no measurable current flowing in the dark and typ- 
ical photocurrents of a few nanoamperes when the excitation 
beam is applied. On this particular device, a larger current 
'leak' existed between the left and lower gates, as shown in 
the lower panel of Fig. S2. When Vlowoi — VLcit ^ 20 V, a 
current of several microamperes flows between them. 

It is important to note, however, that the presence of a 
photocurrent between gates near the diamond surface does 
not preclude the accumulation of a space charge deeper in 
the sample, especially since photoexcited charge transport 
through diamond is an inefficient process, and the timescales 
for charge storage are very long [3]. As illustrated schemati- 
cally in Fig. 2b of the main text, the out-of-plane component 
of the photoinduced field is primarily due to positive charge 
from photoionized donors below the NV center, which results 
from the effective voltage applied across the sample rather 
than between two lateral gates. Thus even with a substantial 
leak between two gates, both the dielectric and photoinduced 
components of the dc Stark shift are always present. Practi- 
cally, however, photocurrents may lead to local heating of the 
sample, and we attribute the leak between the left and lower 
gates to the broadening of the PLE resonances observed in 
Fig. 3 of the main text when larger lateral voltages are applied. 



ANALYSIS OF DC STARK SHIFTS 

The first step in the analysis of dc Stark shifts in our de- 
vices is a determination of the NV center's orientation in the 
diamond crystal. This may be any of the eight projections 
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FIG. S2. Top panel: I-V curves showing conductivity between gate 
the right-hand gate of device A (see Fig. Ic in the main text) and the 
remaining three surface gates, both with (green circles) and without 
(purple squares) 532 nm excitation. Lower panel: Similar photocon- 
ductivity measurements (excitation applied) with the bias connected 
to the left-hand gate (red circles) and lower gate (blue squares), 
showing a larger leakage current between these two gates. 



from the (111) family; note that antiparallel (i.e., N-V com- 
pared to V-N) orientations are distinguishable due to the per- 
manent dipole in the z direction. Using the coordinate systems 
described in the text, we note that N-V axis projects either 
in the ±X or ±Y sample directions, and also has a compo- 
nent along ±Z. We can therefore label a given orientation 
by [px,Py ,Pz), where px,Y = 0,±1 give the sign of the 
N-V projection in the {X, Y) axes (note p\ + Py — 1) and 
Pz = ±1 gives the sign of the Z-projection. For a given NV 
center, we can determine whether px or py equals zero from 
the dependence of the photoluminescence on the linear polar- 
ization of the 532 nm repump beam (the photoluminescence 
is a maximized when the light is polarized perpendicular to 
the N-V axis). It is then usually straightforward to determine 
the exact projection from the slope of the longitudinal dc Stark 
component D in response to applied voltages, keeping in mind 
that the rectified field points mainly out of plane. For exam- 
ple, we know from polarization measurements that the NV 
center of Fig. 1 projects in the zLX direction, and from the 
sign of D in response to Vx (Fig. 2c) we see that px = +1 
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and Pz = —1, corresponding to the crystal direction [111] as tween sample and NV-center coordinate systems is given by 
depicted in the schematic of Fig. 2b. 

Fnv = MF Sample, (SI) 

Given the NV-center symmetry axis, the transformation be- where M is the rotation matrix 



Px Py -V2p2 

M = 4? I -V^PyPz V^PxPz 



V3 



(S2) 



\/2px \/2pY Pz 



Applying this transformation for the NV center in our lateral device, we obtain expressions for the dc Stark components v and 6 
in terms of the local field in sample coordinates: 



1/2 



(S3a) 
(S3b) 



By expressing the local field according to our model as F = 
(iVv + (3\V\^, we obtain a dependence on applied bias that 
we can quantitatively compare to our measurements. 

To reduce the number of free parameters, we fix the 
components of the dielectric unit vector v based on elec- 
trostatic simulations of our device using the COMSOL 
MULTIPHYSICS® software. For a single set of mea- 
surements (e.g., i? and 6 as a function of bias along 
a single direction), this leaves eight free parameters: 
{ujo,SE^,SEy- 0Am,fi±/Am,£,x,£.Y,^z}- This equals 
the combined number of free parameters expected for the form 
of both P and 5 (i.e., two lines and two hyperbolae constrained 
to meet atV = 0) and so results in a well-conditioned mini- 
mization problem. When we have multiple measurements of 
the same NV center, e.g., for the repump/no repump mea- 
surements in Fig. 2d, we can constrain the parameters even 
further, and require that { wq , Se^ , , /?A/X|| , ij,± / A/zy } are 
common for each NV center, while only the rectification vec- 
tor ^ varies independently for each set. In the case of the 
lateral measurements of Fig. 2c, we also allow the dielectric 
coefficient p to vary between the Vx and Vy measurements, 
and the resulting best-fit parameters for a combined fit as de- 
scribed above are listed in Table I. As expected, the rectified 
field points mainly in the direction. 



Electrostatics of the photoinduced fields 

A full description of the photoinduced fields in our devices 
is beyond the scope of this work, especially since our phe- 
nomenological model seems to capture the behavior impor- 
tant for controlling NV-center Stark shifts in these geome- 
tries. Such a description would presumably need to account 
for the gate geometries of our devices, the variation of optical 
intensity through the sample in a confocal geometry, trans- 



port dynamics of photoexcited charges, and potentially spatial 
inhomogeneities in the distribution of deep donors. Further- 
more, although we expect that substitutional nitrogen atoms 
dominate the electrostatics, other defect centers may also play 
an important role, such as vacancy complexes unintentionally 
created during the irradiation and anneaUng process. Based on 
the electron irradiation dose used, we expect a vacancy con- 
centration (before annealing) of x 10^^ cm""^ [4], which 
is roughly an order of magnitude less than the nitrogen con- 
centration in these samples, and these acceptor defects prob- 
ably have an effect on the transport of photoionized charge 
through the sample. Still, we can gather some quantitative 
information about the photoinduced charge distributions from 
the NV-center Stark response, which we can compare to a sim- 
plified 'band-bending' model treating the device as a Schottky 
diode. 

The top-gate geometry of device B is the simplest to ana- 
lyze, and as mentioned in the text we can compare the mag- 
nitudes of the dc Stark shifts observed in that device (Fig. 2a) 
with those of device A (Fig. 2c) to estimate the local electric 
field strength. Assuming that in the top-gate device the local 
electric field points only in the +Z direction, we relate the 
slope of P to the local electric field Fz = PzV through the 
relation 

dP _ /3zA^|| 
dV ~ y/3h ' 



(S4) 



where the factor l/\/3 appears due to the projection of the 
NV-center symmetry axis in the sample. Below the threshold 
bias in Fig. 2a, dP/dV « 0.06 GHz/V. In comparison, from 
the fits to the lateral data in Fig. 2c described in the previ- 
ous section, we find that the dielectric response of device A 
is ,8LateraiA/U||//i w 0.2 GHz/V. Therefore the ratio of the 
dielectric field in device A to the rectified field in device B is 
/^Lateral/ l^z ~ 2, while from dielectric response alone this ra- 
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TABLE I. Best-fit parameter values from fits of our model to the data in Fig. 2c in the main text, for Stark shifts as a function of applied bias 
(Vx, Vy) in orthogonal directions. Uncertainties are at 95% confidence from the fit alone, and do not include systematic errors due to the 
somewhat strong covariance between many of these parameters and the uncertainty in the fixed direction v of the dielectric field. Based on 
electrostatic modeling, we use v = (0.998, 0.04, 0.07) and (-0.002, 0.9998, 0.02) for the Vx and Vy dependence, respectively. 





Common 


Direction-dependent 


Bias direction 




parameters 


parameters 


Vx 


Vy 


LJo /27r 


0.66 ±0.03 GHz 


/3A^i||//i 


0.334 ±0.002 GHz/V 


0.175 ±0.003 GHz/V 


SeJH 


5.06 ±0.03 GHz 


Cx 


0.15 ±0.01 


0.19 ±0.02 


S^Jh 


0.78 ±0.06 GHz 


iv 


-0.43 ± 0.02 


-0.30 ±0.02 


/ix/A/i|| 


1.41 ±0.02 




1.01 ±0.01 


1.60 ±0.02 



tio should be «100, corresponding to an enhancement of the 
vertical electric field below the top gate by a factor of «50. 

If we picture the top-gate device as a Schottky diode, this 
comparison suggests that when negative (reverse) bias is ap- 
plied, the potential drops across a distance of only «10 ^.m 
rather than the full 0.5 mm thickness of the sample. Assum- 
ing uniform ionization of deep donors with density A^^j, the 
equivalent Schottky 'depletion width' obtained from the solu- 
tion to Poisson's equation is given by [5] 




(S5) 



in terms of the Schottky barrier height if and diamond's di- 
electric constant e « 5.7. For Wd ~ 10 |xm and (f w 
10 V, this estimate yields an ionized charge density of A'^ ~ 
6 X 10^'^ cm^'^, which is in agreement with a wlO% ioniza- 
tion of the anticipated substitutional nitrogen concentration of 
a few parts per billion (wl x 10^^ cm~'^). This rough calcula- 
tion shows that the photoinduced fields we observe are consis- 
tent with the expected concentration of deep defects, although 
it does not account for nonlinearities in the distribution of pho- 
toionized charge. For example, we observe a qualitatively 
similar enhancement of the Z-direction field when focussed 
on NV centers up to w50 |J.m below the top gate, suggesting 
that the charge distribution is not confined to a layer near the 
surface gate but rather follows the region of higher optical in- 
tensity. 



The role of the repump cycle on photoinduced fields 

As described in the main text, we performed a control 
experiment to investigate the influence of the 532 nm re- 
pump cycle on the rectified component of the local electric 
field. Using the device shown in Fig. 2d, we apply a mi- 
crowave current through a short-terminated waveguide reso- 
nant with the ground-state \Sz) ^ {\Sx) , \Sy)} transition at 
2.878 GHz. This maintains a mixed spin state during our mea- 
surements since, in the absence of the repump cycle, weak 
spin-nonconserving transitions would otherwise polarize the 
spin away from resonance with the red laser In the PLE spec- 
tra (Fig. S3), we therefore observe the spin-conserving optical 
transitions for all three triplet spin states. We isolate the \Sz) 
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FIG. S3. PLE spectra as a function of gate voltage, for the NV center 
circled in Fig. 2d, both with the standard 532 nm repump cycle (left) 
and with only intermittent repump pulses applied at zero bias (right). 
The \Sz) transition frequencies reconstructed from the fits shown in 
Fig. 2d are overlaid (dashed curves). Other resonances visible in 
these spectra correspond to spin-conserving optical transitions for 
the , 15^)} spin states, as described in the text. 



resonance frequencies from fits to these PLE spectra to com- 
pute the dc Stark components v and 8 plotted in Fig. 2d of 
the main text. These are marked by dashed curves in Fig. S3, 
which are the resonance frequencies reconstructed from the 
best-fit results of our model plotted in Fig. 2d. The analysis 
leading to the fits is similar to that described above (in this 
case the NV-center projection is [111]), and the best-fit pa- 
rameters are listed in Table II. 
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laser on the photoinduced field. As in Table I, imcertainties are quoted at 95% confidence from the fit, but do not account for covariance 
between the fit parameters and the direction of the dielectric field, fixed to v = (0.97, 0.04, —0.24) based on electrostatic modeling. 
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Repump-dependent 








parameters 
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0.06 ±0.01 
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0.30 ±0.03 
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